Spatial Structures in a Generalized Ginzburg-Landau Free Energy 
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Searching for characteristic signatures of a higher order phase transition (specifically of order 
three or four), we have calculated the spatial profiles and the energies of a spatially varying order 
parameter in one dimension. In the case of a p"* order phase transition to a superconducting 
ground state, the free energy density depends on temperature as a'' , where a = ao{l — T/Tc) is 
the reduced temperature. The energy of a domain wall between two degenerate ground states is 
€p ~ a''"^/^. We have also investigated the effects of a supercurrent in a narrow wire. These 



effects are limited by a critical current which has a temperature dependence Jc{T) 



j(2p-l)/2^ 

The phase slip center profiles and their energies are also calculated. Given the suggestion that the 
superconducting transtion in (Bai_a;Ka;)Bi03, for x = 0.4, may be of order four, these predictions 
have relevance for future experiments. 



I. INTRODUCTION 

A large part of statistical mechanics, having to do with 
the physics of phase transitions and nonlinear properties 
is based on a rather fundamental generalization, due to 
Landaui*^, of the well-known Gibbs free energy. There 
is a functional F depending on an (in general complex) 
"order parameter" M(r), given by. 



F[M(r)] = / dv[-a|M|2+6|M|4+c|VMp-M-H], (1) 



where a oc (1 — T/Tc) and h and c are nonnegative, 
temperature independent constants. Here H represents 
an external field which couples directly to the order pa- 
rameter. 

This free energy is the backbone of a vast literature^ 
on second order phase transitions, critical phenomena as 
well as nonlinear properties. Regardless of the micro- 
scopic fabric of the system under consideration, near a 
second order phase transition where universal features 
appear which are independent of the details of interac- 
tion, this free energy contains all of the necessary infor- 
mation. Thus a superconductor or a magnet or a nematic 
liquid crystal, all have the same Ginzburg-Landau (GL) 
free energy as long as the appropriate order parameter M 
is identified. It is however limited to describe only the 
vicinity of a second order phase transition. When the 
fluctuations are ignored, the corresponding Gibbs free 
energy G, identified as the minimum of the functional F , 
has the temperature dependence of G ~ a^. 

According to the classification of phase transitions as 
proposed by Ehrenfesli^, the transitions in general can be 
of any order. In a second order phase transition, the spe- 
cific heat and the compressibility, which are the second 
order derivatives of the free energy (with respect to tem- 
perature and pressure respectively) , are discontinuous at 
the transition. In general though, the discontinuity is 
often replaced by a weak, often logarithmic singularity. 
One could then view the Ehrenfest definition of an order 
as one where the lower derivatives of the free energy are 



continuous at the transition but the higher derivatives 
are singular. Thus in a fourth order phase transition, 
all third order derivatives are continuous, all fifth order 
derivatives are singular at the transition. The fourth or- 
der derivatives are either discontinuous or are weakly sin- 
gular. 

The generalization of the Landau free energy to a 
higher order phase transition is achieved via a weight 
function. For a third and a fourth order phase transi- 



tions, we nee' 
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respectively, within mean field in the 



sense discussed above, G ~ and G ~ a*. This can be 
obtained from "Landau-like" free energies defined with a 
weight function: 

Fill = I dv\M\^[~a^\M\'' + 63|M|4 + calVAfp], (2) 



Fiv = j dv\M\'^[-a4\M\^ + bi\M\'^ + Ci\VM\% (3) 

In general the weight function is given by \M\^^p~'^'^ for 
p > 2. Here a^'s change sign at Tc and are often used 
below for the reduced temperature Op — ap°{l — T/Tc). 
The other constants bp and Cp, p = 3,4 are non- negative 
and generally temperature independent. In a charged 
fluid (such as a superconductori^), the coupling to a 
magnetic field takes the form of a gauge transformation 
V — > V — 2niA/(t)o- Here 0o = h/2e is the superconduct- 
ing flux quantum. The corresponding Gibbs free energy 
is to be identified with the minimum of the functional F 
above with respect to the order parameter. It has the 
expected temperature and field dependence. 

The importance of these free energies is highlighted 
by the recent discoverji^iiS of anomalies in the supercon- 
ducting transition in (Bai_3;K2;)Bi03(x = 0.4). These 
anomalies, a missing discontinuity in specific healiii^ as 
well as in susceptibility'^'^*', indicate a phase transition of 
order higher than second. The actual order in BKBO has 
been determined^ to be fourth. In another example, the 
specific hea^ in Bi2212 (chemical compound)has been 
noted to be kink-like, representing a third order phase 
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transition. The tell-tale signature of a third order phase 
transition, namely a penetration depth ~ remains 
to be investigated. There are other examples in Ref. Q 
where some of the anomalies have been seen but a de- 
tailed study still needs to be carried out. Regardless 
of the detailed microscopic reasons, the free energies in 
Eqs. (|2I3|) characterize the phase transition in its en- 
tirety. For BKBO, the function M is the superconducting 
order parameter. In the absence of a microscopic theory 
we do not know whether it is also the energy gap at the 
Fermi surface, although the tunnelling measurements^^ 
seem to be consistent with this identification. The earlier 
analysi g^i^i'' has been based on a uniform order parame- 
ter, except for some straightforward effects of a magnetic 
field such as a derivation of the London equation and the 
temperature dependence of the penetration length. 

This paper is aimed at studying the properties of a 
substantially inhomogeneous order parameter such as 
found--'ASiiI in defects and textures in the condensed 
state. Their analysis is sometimes carried out within 
homotopy theorjJ^ where the geometrical/topological 
aspects can be studied in detail. These (geometri- 
cal/topological) properties, characteristic of a mapping 
between the real space and the order parameter space, 
will be unaffected by the order of the transition. All 
quantitative properties'^ for example the energies of de- 
fects, how do two defects interact or the effects of the 
defects and textures on the thermodynamic properties; 
they will all be different. In this paper we are interested 
in the solutions of the relevant nonlinear differential equa- 
tion. To begin with, let us limit ourselves to one space 
dimension. The defects discussed here are domain walls 
in three dimensions (3D) or topological point defects that 
occur in narrow wires. In physical terms, these are the 
energies and the profiles of the order parameter between 
two degenerate superconducting states. A special case 
in Sec. IV below refers to phase slip centers. Voltage 
drops appear across these objects in a narrow wire in the 
presence of a current. 

To recall the results of an earlier analysis^iSii for a p*'* 
order phase transition: (1) The thermodynamic free en- 
ergy at the order parameter minimum follows Fop oc a^. 
(2) The superfluid density satisfies ps oc a^"^. This 
leads to the result that ~ Hd oc a^"^, where A 
is the London penetration depth. (3) Since the coher- 
ence length ^, which measures the stiffness of the order 
parameter, has the temperature dependence oc a, the 
Landau parameter k — X/^ is temperature dependent (in 
contrast to a superconductor undergoing a second order 
phase transition where «; is a temperature independent 
constant) and follows, k oc a^~P/'^. 

The principal results in this paper are: When a defect 
is created, the order parameter is suppressed in a small 
region of order of the coherence length, ^. If the bulk 
condensation energy density (energy per unit length) is 
Eo, the energy of a defect is oc E'oC- This has the tem- 
perature dependence of a^"^^^. The numerical factor 
in front of this expression is calculated below in detail 



for some specific cases. The critical current in a nar- 
row wire is known to be JciT) oc a' with g = 3/2 for a 
second order phase transition. Below we derive this ex- 
ponent for a third and a fourth order phase transition, 
where respectively q — 5/2 and 7/2. For a pih. order 
phase transition, the exponent can be estimated as fol- 
lows. Consider the free energy of a narrow wire in the 
presence of a current J; it must be F ~ ~ /2ps- It 
follows readily that given the temperature dependences 
described above, q = {2p — l)/2. The energy of a phase 
slip center, and therefore the temperature dependence 
of resistivity near a suprconducting phase transition was 
first calculated by Langer and AmbegaokariSii^. A ther- 
mal distribution of these localized voltage points leads 
to an activated/exponential rise in resistance with the 
activation energy being the energy of a single phase slip 
center. The latter in turn is the above mentioned energy 
of a defect, of the order of oc Eo£,. 

The outline of this paper is as follows. Section II con- 
tains the mathematical formalism we have used. In par- 
ticular this section contains the Euler-Lagrange equation 
for a one dimensional nonlinear field theory derived from 
Eqs. H2I3|I . In Sec. Ill, we describe the solutions of 
the nonlinear partial differential equations. The current 
induced effects are discussed in Sec. IV. These include 
temperature dependence of the critical current in a nar- 
row wire and the temperature dependence of the resis- 
tance near the phase transition. Finally the last section 
contains a summary of our conclusions. 



II. MATHEMATICAL INFRASTRUCTURE 

The gradient terms in Eqs. (|2I3|I contain information 
about the stiffness of the order parameter. The precise 
degree of this stiffness, represented by the coefficients 
and which can be measured in a superconductor in ei- 
ther the upper critical field (the resistance of the order 
parameter to spatial variation) or in the London pene- 
tration depth (and the lower critical field, measuring the 
resistance to magnetic field), depends on the material 
parameters. The limitation to the lowest order gradient 
terms is initially motivated by an esthetic curiosity about 
the long wavelength phenomena. Eventually though, the 
ability of this formalism to describe defects and textures 
is an a posteriori rationalization. 

An Euler-Lagrange equation for the free energy (for a 
scalar order parameter M) described in Eq. (0) is given 
by 

= -2a3M^ + 3fe3M^-C3[M(VM)2-|-M2(V^M)]. (4) 

In one dimension (ID), V = i9z and = dzz- The 
letter subscripts denote differentiation. The dimensional 
variables can be scaled out with the transformations 

^-^^\/^"! 2/ = 2/6; 6^ = 03/203; (5) 
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Eq. Q then becomes 

The function f ^ v? satisfies 



0. 



fyy -f + f^O. 



(6) 



(7) 



The corresponding energies are then calculated to be 

(2a3)' 



Fill — £"363; i?3 — £,3 



(353)2 ' 



and 



£3 
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dy 



3 fdf 
2 



1 - if + 2f 



(8) 



(9) 



Since the defects are localized over a spatial region of 
order ^3, we notice that the energy scale follows the tem- 
perature dependence 03^3. The function n(y) approaches 
±1 in the bulk. The dimensionless energy expression has 
been adjusted to define the defect energy with respect to 
the state with a uniform order parameter. 

The functional has to be limited to f{y) > 0. 

If it were not so, then Eq.@ as a functional of / would be 
unbounded from below and the mean field theory (where 
the sum over possible configurations is believed domi- 
nated by the minimum energy configuration) is then un- 
defined. Since f ^ v?, the condition is also a natural 
requirement. 

The physics derived from Eq. ((3J) for a fourth order 
transition is similar, although analytically less tractable. 
The Euler-Lagrange equation, in this case, becomes, 

Ci[M^M^^ + 2M^{M^f] + ia^M^ - Ah^M'^ = 0. (10) 

Following very similar transformations (removal of di- 
mensional variables), but with an important difference, 
namely that / cx M^, the equation corresponding to Eq. 
O becomes; 



\fvv + .f- = 0, 



(11) 



where f/^ = n ^ ^^/y tft and y = z/^4 with = 

£4/304. Also, the energy of these textures is given by 
Fiv = Ei€i with = ^4(3a4)'*/(464)3. The quantity 
£4, as an integral, is given by 

1 - 4/2 + 3/^/3 + 



£4 = 



dy 



dl 
dy 



(12) 



III. SPATIAL CONFIGURATIONS 

In this section we calculate the spatial profiles. This 
is done separately for the third and fourth order free en- 
ergies. The solutions are for a scalar order parameter in 
one dimension. 



A. Third Order Free Energy 

In case of Eqs. (0 and the transformation f = n? 
eliminates a sign degeneracy. It also restricts / > 0. 
A domain wall between regions with n = 1 and n = 
— 1, the two degenerate ground states in bulk, becomes 
a profile where f — 1 everywhere, except in a narrow 
region of order ^, where it vanishes. Let us first consider 
the solution for Eq. {T)). A first integral is obtained by 
multiplying the partial differential equation by ny and 
integrating over y. This leads to the first integral, 
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K. 



(13) 



Recall that the traditional analysis from here onwards, 
notices, that if y were time, K would be the energy of a 
particle with position /, moving in a potential well —V, 
where V is the gradient-free part of the free energy (for 
/ > 0). For small K, the solutions are periodic, similar to 
a periodic arrangement of holes in the condensate. Equa- 
tion (O may be integrated in terms of elliptic functions. 
As K increases, the sinusoidal functions sharpen into 
square wave like structures, their wavelength increases 
until eventually for K Kq — 1/6, one obtains a soli- 
tary wave solution. (In general, Kq — l/2p.) The center 
of the solitary wave in f{y) is linear. Since at the center 
the order parameter is small it follows that f{y) ~ 2y/l{y. 

For small K, the potential for the motion of the fic- 
titious particle is simply f. The motion in a quadratic 
well is described by / = V2K\ sin(?/\/2)|. Since / > 0, 
the analytic solution needs to be understood in light of 
this positivity constraint. 

The energy of the small amplitude periodic structure 
(up to order linear in K) is zero, as can be evaluated by 
direct calculation. As K increases and the periodic struc- 
ture begins to resemble a soliton anti-soliton lattice with 
wide separation between the kinks, the energy reduces to 
the sum of the rest masses of the individual kinks, re- 
duced by the attractive interaction between them. For 
large distances d between two domain walls, the asymp- 
totic interaction between a soliton and an antisolitoniS 
is given by 



U{d) = -36\/2 e'^"^. 
Integrating Eq. (|13|1 for AT = 1/6 we get 



/ 



3 tanh 



X 



V2 



(14) 



(15) 



Here X is the second integration constant, it represents 
the translational freedom in locating the center of the 
solitary wave anywhere. Since / > 0, it follows that 
X > = \/2tanh"^ ^ = 0.93. The solution for y < 
is obtained by folding the result around the vertical axis. 
Figure. (J^l shows the solution for f{y). This is different 
from a solution one obtains for a conventional solitary 



FIG. 1: The function f{y) as a function of y for a third order 
free energy. 



FIG. 3: The approximate solution, Eq. 1)20^ (sohd hne), and 
numerically exact solution (crosses) of the function f{y) = 
for a fourth order phase transition, obtained from Eg. 1)181 



wave as in Eq. for a second order transition. There 
is a discontinuity of slope at y = which is remnant of 
the similar discontinuity in the small amplitude solutions. 
It arises from the constraint / > 0. 

The solution of the original Eq.lO is obtained from 
n — ^/f with the positive root y > and the negative 
root for y < 0. The function n{y is shown in Fig.Q. 

The energy of the solitary wave structure is given by 



£3 



3^2 



3\/3 



= 0.522. 



(16) 



The n(?/) solution is depicted in Fig. |5J). This is gen- 
erated from the solution for f{y) shown in Fig. ^ by 
taking the positive square root for y > and the nega- 
tive square root for y < 0. Near the center the function 




FIG. 2: Domain wall profile n{y) for a third order phase tran- 
sition. 



B. Fourth Order Free Energy 

For a fourth order transition free energy, Eq. the 
Euler-Lagrange equation, Eq. Hll(l . can be readily inte- 
grated once. The resulting expression is 
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ify 



-J - of 



8/3 



K. 



(17) 



Again for small if, one has the periodic solutions which 
start out as sinusoidal solutions but sharpen into square 
wave like functions as K increases. For K ~ 1/8, the 
solution for a single defect is given implicitly by the in- 
tegral 



df 



V3 



(18) 



The solution /(y) with /(O) = 0, is sketched in Fig.© 
as the set of crosses. A numerical integration of Eq. 1)12(1 
leads to the value of the integral as 64 — 0.241. The defect 
energies get smaller for the higher order free energies. 
The ratio 63/64 = 2.17. In comparison, the domain wall 
energy for Eq.©, 63 = 2\/2/3 = 0.943 and 62/63 = 1.81. 

An approximate solution for Eqs. pOI lll|l can be 
worked out by starting with the free energy in dimension- 
less variables, as described in Ea. (ll2ll . This free energy 
can be broken into two parts; t^^ — A + B where. 



A = 



12 



dy 



1 - 3/2 + 2f 



dl 
3 \dy 



and 



B-^J dy[3f^^-f{l + 2f)]. 



(19) 



(20) 
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The separation here is guided by the solvabihty. Its va- 
hdity is justified a posteriori. The unperturbed part of 
the free energy, represented by A is similar to Eq. 0. 
There is an overaU factor of 2 (£4 is smaller) and the size 
of the structure is slightly smaller by a factor -\/8/9. The 
perturbation B is non-zero only when / 7^ 1. 
The solution to Eq. (|19|l is readily obtained. 



r(y) = i[3tanh2(3(j/ + X)/4)-l]. (21) 

The energy of this structure is given by ^ = 0.246, rather 
close to the exact value for the energy (0.241) cited fol- 
lowing Eq. H18|) . The integral for B can be evaluated 
using f°{y). It is equal to —0.001. Shown in Fig.Q as 
the solid line, is the solution described in Ea. (|21|l . The 
difference between the exact numerically integrated solu- 
tion and the approximate one is negligible. 



IV. CURRENT FLOW EFFECTS 

Consider a narrow wire (lateral dimensions smaller 
than the coherence length) made of a material described 



by the generalized free energies. If the wire carries cur- 
rent J, the superfluid density is reduced, first quadrati- 
cally but as the current reaches the critical current, more 
precipitiously. The superfluid density vanishes at a crit- 
ical current Jc(T). However even for J < Jc(T), voltage 
drops appear in small regions called phase slip centers. 
A comprehensive and lucid discussion of both the math- 
ematical as well as physical effects has been provided by 
Langer and Ambegaokar— . 



The physical effects remain qualitatively the same for a 
higher order phase transition. The mathematical details, 
however, are different for the free energies in Eqs. 

In the presence of a constant current J, the free en- 
ergy corresponding to Eq. ||2Jl for a third order transition 
becomes. 



J 



F,^F- J.A 



dvlMl 



-asm 



b.\M\ 



4C3M6 



(22) 



Here A is a. vector potential, conjugate to the super- 
current J. Because of the use here of the Legendre's 
transformation to develop Fj as a function of specified 
J, the final equation is entirely in terms of the amplitude 
of the order parameter. Its minimum yields the equi- 
librium order parameter, including the effects of a finite 
supercurrent J. 

The physics described in Eq. (|22|l begins with the usual 
free energy in Eq. for J = 0. With increasing J, the 
order parameter M corresponding to the minimum of the 
free energy moves to smaller values. The free energy also 
has a maximum which separates the minimum from the 
normal state for M = 0. Both the minimum as well as 
the maximum disappear for J > Jc{T) and there is no or- 
der parameter solution. The minimum free energy order 
parameter solution M{J) is a solution of the equation. 



j2 = c3M^[4a3M^ - 663 Af^]. 



(23) 



For each J < Jc, there are two solutions. The smaller 
M corresponds to the local maximum of the free energy. 
The larger solution closer to the zero current equilib- 
rium value is the global minimum and ceases to exist 
for J > Jc{T). The corresponding problem for Eq. 
leads to a temperature dependence of Jc{T) — c?!"^ . For 



Eq. @, the result is Jc(T) ~ 03^/2 _ ^oy Eq. ©, the 
critical current has an even weaker temperature depen- 
dence, Jc(T) ~ 04''/^. It is straightforward to extend 
this calculation to an arbitrary order phase transition. 
For a pth order transition the exponent for the temper- 
ature dependence is (2p — l)/2. We thus have another 
characteristic signature of the order of a transition in the 
temperature dependence of the critical current in a wire. 

The solution corresponding to a phase slip center for 
Eq. was originally obtainedi^ by Langer and Ambe- 
gaokar. In dimensionless variables, Eq. H22I) transforms 
into. 



Fill — E3 



dy 



2/ +3/ +4/^ 2/2 



(24) 



The notation is similar to that following Eqs. (jSJ and 
© with the addition for the current. Thus — 
( J^/2c3)((36)^/(2a3)^). In terms of the dimensionless 
variables, Eq. becomes 



5' = /*(!-/)■ 



(25) 
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The Euler-Lagrange equation is given by 

1 2 



r,fyv 



(26) 



For a given current J (less than Jc{T)), there are two 
solutions to Eq. H25|l . The large f corresponds to the 
order parameter reduced in the presence of the current. 
The smaller f solution corresponds to the maximum in 
the free energy. Equation can be integrated once to 
lead to: 



:fy^ + U{f)^K- U{f) = 



P 



3 



2/2- 



(27) 



Here too, the solutions are those corresponding to the 
motion of a particle with position /, time y, moving in 
a potential U{f). The allowed values of K are bounded 
by Ki < K < K2- Here Ki corresponds to the value 
of U at the local minimum. In the free energy this is 
the maximum corresponding to the smaller / solution of 
Eq. 12511 . For K > Ki, the solutions oscillate about 
the extremum fmax with a wavelength determined by 
other parameters of the problem. K2 referes to the value 
of U at the local maximum. This corresponds to the 
free energy minimum and the order parameter as reduced 
by the current. For K ^ K2, the order parameter / is 
constant everywhere except in a small region, where it 
drops down to a value fmin such that U{fmin) — K2- 

Equation H27I) needs to be finally integrated. The so- 
lution can be expressed as an integral. 



/: 



fdf 



^2p - 3/4 + &Kp - 3.92 



(28) 



The profile is shown in Fig.Q). Qualitatively the pro- 
file is similaoi^ to the one found for Eq. The details 



0.4 ■ 

0.3 ■ 



g=0.035 
g=0.134 
9=0.250 



FIG. 4: Order parameter profile fsiy) for a third order phase 
transition in the presence of a current with g = 0.035, 0.134 
and 0.250. 




FIG. 5: The energy of a phase slip center, numercailly in- 
terated from Eg. 12411 for a third order phase transition free 
energy. 



however are all quite different. Figure ||SJ) shows the en- 
ergy of the phase slip centers as a function of the dimen- 
sionless supercurrent g for the third order free energy. 
When properly scaled, the results for a fourth order free 
energy look identical. 

A similar analysis can also be carried out for Eq. ||2J), 
generalized to include the effect of a supercurrent. In the 
dimcnsionlcss units the free energy becomes. 



IV 



3/ +4/ +9/!. 3J2 



(29) 



Here the notation is identical to that following Eq. lll|l . 
The dimensionless current = (3 j2/4c4)((46)6/(3a)^) 
is given by, 



(30) 



One straightforward consequence of this equation (or 
Eq.(j2ni)) is that the suppression of the order parameter 
due to a supercurrent is given here by p''^ ~ This 
and the result corresponding to Eq. (^5)1 / ~ 1 — look 
different only because / depends differently on the order 
parameter. Here the maximum in g^ occurs (when the 
two solutions merge and there is no solution for larger g) 
at fo = (6/7)3/2 ^ 0.78, corresponding to 5^ = 0.24. 

We have not integrated the Euler-Lagrange equation 
for a phase slip center in this case. In earlier calculations 
related to Eq. (Q, the final calculation seems to be pos- 
sible only numerically. Qualitatively the main features 
of the solutions for Eq. ((21 remain useful in offering in- 
sight into the solutions here. We expect the profile of a 
phase slip center to follow the basic features of the so- 
lution described in Eq. H28() . The phase slip center will 
again appear as a spatially localized suppression of the 
order parameter, over a length scale which will vary with 
current. Overall the energy of this structure will be E^^. 
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V. SUMMARY AND CONCLUSIONS 

In view of the success of a Landau free energy in de- 
scribing issues connected with a second order phase tran- 
sition, we have been stimulated to explore the properties 
of defects and textures in the presence of a higher or- 
der phase transition using generalized Landau free energy 
functionals. We expect that, regardless of the so far un- 
known microscopic details, properties near a higher order 
transition can be described by a generalized free energy 
in essentially the same way that the usual Landau free 
energy captures the essential properties near a second 
order phase transition. 

For a third order free energy, the transformation 
f{y) = n'^iu) changes Eq.Q into Eq.®. In the lat- 
ter version, the order parameter must be positive also 
because otherwise the free energy is unbounded below. 
That would eliminate the notion of exp(— /) as a config- 
uration probability. The constraint however leads to a 
very different mathematical structure for the defect pro- 
file. In particular, the profile has a discontinuity in the 
slope fy at the origin. It is not altogether surprising that 
the positivity constraint leads to the discontinuity. 

For Eq.Q, the transformed free energy is seen in 
Eq.(22Il- The term /^/^ could be defined so that one 



is always taking the positive root. There are, it seems 
no symmetry considerations, which will restrict / and in 
this case one is left to speculate freely. The microscopic 
theories will have more to say about whether the pow- 
ers of the terms in the free energy are justifiable and the 
protocol for the proper approach to / = 0. The con- 
straint / > has consequences for stability of the spatial 
textures which will be reportec^ separately. 

There are defect and current related characteristic sig- 
natures of a higher order phase transition. The current 
induced suppression of the order parameter is quadratic 
(in supercurrent) for small currents. The critical cur- 
rent, at which the order parameter vanishes and super- 
conductivity disappears, has a temperature dependence 
Jc(r) ~ a^^^"-^)/^. We have also calculated the energy 
and spatial profiles of structures such as an interface be- 
tween degenerate superconducting states and a phase slip 
center which appears when the supercurrent approaches 
the critical current. 
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